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1. Introduction 

Knowing whether the orbits of a dynamical system are ordered or chaotic is fundamental 
for the understanding of the behavior of the system. In the dissipative case, this 
distinction is easily made as both types of motion are attracting. In conservative 
systems, however, distinguishing between order and chaos is often a delicate issue (e.g. 
when the chaotic or ordered regions are small) especially in systems with many degrees 
of freedom where one cannot easily visualize the dynamics. For this reason it is of 
great importance to have quantities that determine if an orbit is ordered or chaotic, 
independent of the dimension of its phase space. 

The well-known and commonly used method for this purpose is the evaluation 
of the maximal Lyapunov Characteristic Exponent (LCE) a%. If o~\ > the orbit is 
chaotic. Benettin et al. pQ studied theoretically the problem of the computation of all 
LCEs and proposed in [2] an algorithm for their numerical computation. In particular, 
o"i is computed as the limit for t — > oo of the quantity 

L t = ]■ In , i.e. a x = lim L t , (1) 

t \w(0)\ t-*oo 

where w(0), w(t) are deviation vectors from a given orbit, at times t = and t > 
respectively. The time evolution of w is given by solving the so-called variational 
equations (see Sec. EJ). Generally, for almost all choices of initial deviations u>(0), 
the limit for t — > oo of Eq. (Q) gives the same o~\. 

In practice, of course, since the exponential growth of wit) occurs for short time 
intervals, one stops the evolution of wit) after some time 7\, records the computed 

normalize vector w(t) and repeats the calculation for another time interval T 2 , etc. 
obtaining finally o\ as an average over many Tj, i = 1, 2, . . . , iV as 

N 



N 1 



ai N 

%=i 

The basic problem of the computation of o\ is that, after every Tj, the calculation starts 
from the beginning and may yield an altogether different than the T(i-i) interval. 
Thus, since 0\ is influenced by the whole evolution of w(0), the time needed for L t (or 
the Lrp.) to converge is not known a priori and may become extremely long. This makes 
it often difficult to tell whether o~\ finally tends to a positive value (chaos) or converges 
to zero (order). 

In recent years, several methods have been introduced which try to avoid this 
problem by studying the evolution of deviation vectors, some of which are briefly 
discussed in Sec. EI In the present paper, we focus our attention on the method of 
the Smaller Alignment Index (SALI) 3J, performing a systematic study of its behavior 
in the case of autonomous Hamiltonian systems with 2 (2D) and 3 (3D) degrees of 
freedom. This method has been applied successfully to several 2-dimensional (2d) 
and multidimensional maps [3], where SALI was found to converge rapidly to zero for 
chaotic orbits, while it exhibits small fluctuations around non-zero values for ordered 
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orbits. It is exactly this "opposite" behavior of the SALI which makes it an ideal 
indicator of chaoticity: Unlike the maximal LCE, it does not start at every step a new 
calculation of the deviation vectors, but takes into account information about their 
convergence on the unstable manifold from all the previous steps. The method has 
already been used successfully as a chaos detection tool in some specific Hamiltonian 
systems |U El El El El , although some authors [313 use different names for the SALI. 

The paper is organized as follows: In Sec. E]we recall the definition of the SALI and 
present results distinguishing between ordered and chaotic motion in 2 and 3-degrees of 
freedom (2D and 3D) Hamiltonians, comparing also the efficiency of the SALI with the 
computation of o~\. In Sec. 0]we explain the behavior of the SALI for ordered and chaotic 
orbits, showing that in the latter case SALI converges exponentially to zero following 
a rate which depends on the difference of the two largest Lyapunov exponents o~\ and 
o 2 . In Sec. 0]we demonstrate the ability of the method to reveal the detailed structure 
of the dynamics in the phase space. In Sec. El we compare the SALI method with some 
other known methods of chaos detection and in Sec. El we summarize our results. 



2. Application of the SALI in Hamiltonian systems 



The basic idea behind the success of the SALI method [3 is the introduction of a simple 
quantity that clearly indicates if a deviation vector is aligned with the direction of 
the eigenvector which corresponds to the maximal LCE. In general, any two randomly 
chosen initial deviation vectors u>i(0), ?Z>2(0) will become aligned with the most unstable 
direction and the angle between them will rapidly tend to zero [2j. Thus, we check 
if the two vectors have the same direction in phase space, which is equivalent to the 
computation of the above-mentioned angle. 

More specifically, we follow simultaneously the time evolution of an orbit with initial 
condition x(0) and two deviation vectors with initial conditions Wi(0), w 2 (0). As we 
are only interested in the directions of these two vectors we normalize them, at every 
time step, keeping their norm equal to 1. This controls the exponential increase of the 
norm of the vectors and avoids overflow problems. Since, in the case of chaotic orbits 
the normalized vectors point to the same direction and become equal or opposite in 
sign, the minimum of the norms of their sum (antiparallel alignment index) or difference 
(parallel alignment index) tends to zero. So the SALI is defined as: 



SALI(t) = min 



wi(t) w 2 (t) 



+ 



Wl(t) w 2 (t) 



\Wi 



\w 2 



(2) 



\Wl(T)\\ \\W 2 ( 

where t is the time and || - || denotes the Euclidean norm. From the above definition is 
evident that SALI(t) e [0, y/2] and when SALI = the two normalized vectors have the 
same direction, being equal or opposite. 

In order to apply the SALI method to Hamiltonian systems, we shall use here two 
simple examples with 2 and 3 degrees of freedom: the well-known 2D Henon-Heiles 
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system ^Uj, having the Hamiltonian function: 

H2 = \{pl + p\) + \{x 2 + y 2 ) + x 2 y - ^y 3 , (3) 
with equations of motion: 

x = —x — 2xy , y = —y — x 2 + y 2 , (4) 
and the 3D Hamiltonian system: 

#3 = \{vl + P 2 y + Pi) + \{Ax 2 + By 2 + Cz 2 ) - exz 2 - rjyz 2 , (5) 
with equations 

'x = —Ax + tz 2 , y = —By + r/z 2 , z = —Cx + 2z(ex + r]y) , (6) 

studied in ^T] ^] . We keep the parameters of the two systems fixed at the energies 
H 2 = 0.125 and H 3 = 0.00765, with A = 0.9, B = 0.4, C = 0.225, e = 0.56 and 77 = 0.2. 

A simple qualitative way of studying the dynamics of a Hamiltonian system is by 
plotting the successive intersections of the orbits with a Poincare surface of section (PSS) 
. This method has been extensively applied to 2D Hamiltonians, as in these systems 
the PSS is a 2-dimensional plane. In 3D systems however, the PSS is 4-dimensional and 
the behavior of the orbits cannot be easily visualized. One way to overcome this problem 
is to project the PSS to spaces with lower dimensions (see e.g. [T4"l I15j). However, even 
these projections are often very complicated and difficult to interpret. 

In order to illustrate the behavior of the SALI in 2D and 3D systems we first consider 
some representative ordered and chaotic orbits. In Fig. ^a) we plot the intersection 
points of an ordered and a chaotic orbit of Eqs. (J1J), with a PSS defined by x = 0. The 
points of the ordered orbit lie on a torus and form a smooth closed curve on the PSS. 
On the other hand, the points of the chaotic orbit appear randomly scattered. The time 
evolution of the SALI for these two orbits is plotted in Fig. Hfb). In the case of the 
ordered orbit (solid line) the SALI remains different from zero, while in the case of the 
chaotic orbit (dashed line), after a small transient time, the SALI falls abruptly to zero. 
At t ~ 800 the SALI becomes zero as it has reached the limit of the accuracy of the 
computer (10 -16 ), which means that the two deviation vectors have the same direction. 
Thus, after t pa 800 the two normalized vectors are represented by exactly the same 
numbers in the computer and we can safely argue, that to this accuracy the orbit is 
chaotic. Actually, we could conclude that the orbit is chaotic even sooner, considering 
that the directions of the two vectors practically coincide when the SALI reaches a small 
value, e.g. 10~ 8 after some 400 time units. Entirely analogous behavior of the SALI 
distinguishes between ordered and chaotic orbits of the 3D Hamiltonian (J3J), as shown 
in Fig. El 

The initial deviation vectors, w = (dx,dy,dp x ,dp y ), used for both orbits of Fig. 
[U are u>i(0) = (1,0,0,0), 102(0) = (0,0,1,0), but in general any other initial choice 
leads to similar behavior of the SALI. The validity of the above statement is supported 
by the following computations. Focusing our attention on the more interesting case of 
chaotic motion we study the chaotic orbit of Fig. ^ by fixing one of the initial deviation 
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Figure 1. (a) The PSS of an ordered and a chaotic orbit with initial conditions x = 0, 
y = 0.1, Vx ^ 0.49058, p y = and x = 0, y = -0.25, p K ~ 0.42081, p y = respectively, 
for the Henon-Heiles system ©. The ordered orbit corresponds to a closed (solid) 
elliptic curve, while the chaotic one is represented by the dots scattered over the PSS. 
(b) The time evolution of the SALI for the two orbits of panel (a) in log-log scale. The 
solid line corresponds to the ordered orbit while the dashed line corresponds to the 
chaotic orbit. 




Figure 2. The time evolution of the SALI for an ordered (solid line) and a chaotic 
orbit (dashed line) of the 3D Hamiltonian (J5J, with initial conditions x — 0.01054, 
y = 0.05060, z = 0, p x = 0, p y = 0, p z ~ 0.11906 and x = -0.07310, y = 0, z = 0, 
p x = 0.07695, p y — 0, p z ~ 0.06760 respectively. 
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Figure 3. The time T needed for the SALI to become less than 1CP 12 in the case of 
the chaotic orbit of Fig. ^ when we use as initial deviation vectors wi(0) = (0,1, 0, 0) 
and (a) w?2(0) = (0,cos#, 0,sin#), and (b) W2(0) being a random vector. T is plotted 
as a function of 9 in (a), and as a function of a counter i of the randomly generated 
vectors in (b). 

vectors, keeping it e. g. Wi(0) = (0, 1, 0, 0) and varying the second one 102(0). For every 
different pair of initial vectors w\(0), u>2(0) we compute the time T needed for the SALI 
to become smaller than a very small value e. g. 10~ 12 and check if the value of T depends 
on the particular choice of initial deviation vectors. We choose 102(0) in two different 
ways. Firstly we consider the two initial vectors to be on the PSS of Fig. having an 
angle 9 between them so that to 2 (0) = (0, cos#, 0, sin 9). In Fig. 0(a) we plot T as a 
function of 9 for 9 G [0,7r]. As expected T = for 9 = and 9 = it since the two 
vectors are initially aligned. The maximum value of T for 9 ~ 0.4057T corresponds to 
the case of 102(0) being almost perpendicular to the unstable manifold which passes 
near the initial condition of the orbit. In this case the component of 10*2(0) along the 
unstable direction is almost zero and thus, the time needed for the vector to develop a 
significant component along this direction, which will eventually lead it to align with 
the other deviation vector, is maximized. We remark that for all 9 e (0, n) T does not 
change significantly, as it practically varies between 400 and 500 time units. As there is 
no reason for wi(0), w 2 (0) to be on the PSS the second test we perform is to compute 
T for 1000 ^2(0) whose coordinates are randomly generated numbers (Fig. 03). From 
the results of Fig. El we see that T practically does not depend on choice of the initial 
deviation vectors. 

On the other hand, the computation of the maximal LCE, using Eq. (|Tj) . despite 
its usefulness in many cases, does not have the same convergence properties over the 
same time interval. This becomes evident in Fig. 0] where we plot the evolution of the 
L t (panel (a)) and the SALI (panel (b)) for a chaotic orbit of Eqs. (p£J). At t « 1900 the 
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Figure 4. The time evolution of (a) the L t and (b) the SALI for the chaotic orbit 
x = 0, y = -0.01597, p x ~ 0.49974, p y = of the 2D system ©. 

SALI reaches the value 10 -16 and no further computations are needed. Of course we 
could be sure for the chaotic nature of the orbit before that time, for example at t « 1000 
(where the SALI « 10~ 8 ). On the other hand, the computation of the L t (Fig. EJa)) 
up to t ~ 1000 or even up to t ~ 1900, still shows no clear evidence of convergence. 
Although Fig. Ufa) suggests that the orbit might probably be chaotic, it does not allows 
us to conclude its chaotic nature with certainty, and so further computations of L t are 
needed. Thus, it becomes evident that an advantage of the SALI, with respect to the 
computation of L t , is that the current value of the SALI is sufficient to determine the 
chaotic nature of an orbit, in contrast to the maximal LCE where the whole evolution 
of the deviation vector affects the computed value of L t . 

3. The behavior of the SALI for ordered and chaotic motion 

As we have seen for ordered orbits the SALI remains different from zero fluctuating 
around some non-zero value. The behavior of the SALI for ordered motion was studied 
and explained in detail in the case of a completely integrable 2D Hamiltonian ^Hj, in 
which no chaotic orbits exist. It was shown that any pair of arbitrary deviation vectors 
tend to the tangent space of the torus, on which the motion is governed by 2 independent 
vector fields, corresponding to the 2 integrals of motion. Thus, since W\(t) and ^(t), in 
general have one component 'along' and one 'across' the torus, there is no reason why 
they should become aligned and thus typically end up oscillating about two different 
directions. This explains why the SALI does not go to zero in the case of ordered motion. 

Now let us investigate the dynamics in the vicinity of chaotic orbits of a 
Hamiltonian system of n degrees of freedom. An orbit of this system is defined by 
x = (<?i, <?2, • ■ • , Qn, Pi, P2, Pn), with q h p h % = l,...,n being the generalized 
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coordinates and the conjugate momenta respectively. The time evolution of this orbit 
is given by Hamilton's equations of motion 

dx fdH dH\ 



dt \ dp dq 

Solving the variational equations about a solution of (JJJ), x(t), which represents our 
reference orbit under investigation, 

^§ = M(x (t))w, (8) 

where M = ^ is the Jacobian matrix of V, we get the time evolution of an initial 
deviation vector w(to), for sufficiently small intervals [to, t + At}. Note now that, in this 
context, the eigenvalues Ai ^ A2 ^ • • • ^ A 2n of M, at t — t Q , may be thought of as local 
Lyapunov exponents, with e\, e 2 , . . . , e 2n the corresponding unitary eigenvectors. These 
eigenvalues in fact oscillate about their time averaged values, o~\ ^ a 2 ^ ■ ■ ■ ^ o~2 n , which 
are the global LCEs of the dynamics in that region. As is well-known in Hamiltonian 
systems, the Lyapunov exponents of chaotic orbits are real and are grouped in pairs of 
opposite sign with at least two of them being equal to zero [T3]. Thus, the evolution of 
any initial deviation vector W\(0) is given by: 

2n 

^l(t)=5>* WeA< ^> ( 9 ) 
i=l 

where the c*p are in general complex numbers and Aj, depend on the specific location 
in phase space, x(t ), through which our orbit passes. Note that we consider here only 
real eigenvalues and hence real c\ . If some of the Aj are complex, their corresponding 
contribution to will be oscillatory and will not affect the argument that follows. 

To get a first rough idea of the way the SALI evolves, we now make some 
approximations on the evolution of this deviation vector. First let us assume that the Aj 
do not fluctuate significantly about their averaged values and hence can be approximated 
by them, i.e. \ « a^. Secondly, we consider that the major contribution to W\(t) comes 
from the two largest terms of Eq. ©, so that: 

u?i(t) « cS X) e Al * h + 4 1} e Aai e 2 « e ait e x + c? e ff2 ' e 2 . (10) 

In this approximation, we now use (fTU|) to derive a leading order estimate of the ratio 

Ml ~ c?e^e 1+ c?e^e 2 _ cf_ {<n - n)t ( . 

and an entirely analogous expression for a second deviation vector: 

Ml ~ cj^h +cfe^e 2 _ cf { ^ a2)t , , 

w^t - S2ei + o^ e e2 ' (12) 



where Sj = sign(ci ), i — 1, 2. 
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Figure 5. (a) The evolution of L t for the chaotic orbit with initial condition x = 0, 
y = -0.25, p x ~ 0.42081, p y = of the 2D system ©. (b) The SALI of the same orbit 
(solid line) and the function e~ ait (dashed line) for cr\ = 0.047. Note that the t-axis 
is linear. 



In order to compute the SALI, as denned by (j2J) we add and subtract Eqs. fTTf 
and (fl"2"j) keeping the norm of the minimum of the two evaluated quantities. Thus, e\ 
does not appear in the expression of the SALI, which becomes 



SALI(t) = min 



W\(t) w 2 (t) 



± 



± 



c (2) 



.(2) I 



-(CTI— <T 2 )t I 



e 2 || (13) 



Wi(t)\\ \\ir 2 [ 

Denoting by c the positive quantity on the r.h.s of the above equation and using the 
fact that e 2 is a unitary vector we get 

SALI(t) « ce-*" 1-175 '' . (14) 

Eq. (fTlj) clearly suggests that the SALI for chaotic orbits tends to zero exponentially 
and the rate of this decrease is related to the two largest LCEs of the dynamics. 

Let us test the validity of this result by recalling that 2D Hamiltonian systems have 
only one positive LCE o~i, since the second largest is a 2 = 0. So, Eq. (jTlj) becomes for 
such systems: 

SALI(t) ^ce- CT1 *. (15) 

In Fig. 02a) we plot in log-log scale L t as a function of time t for a chaotic orbit of 
the 2D system ©. L t remains different from zero, which implies the chaotic nature of 
the orbit. Following its evolution for a sufficiently long time interval to obtain reliable 
estimates (t ~ 10000) we obtain a x w 0.047. In Fig. \Mb) we plot the SALI for the same 
orbit (solid line) using linear scale for the time t. Again we conclude that the orbit is 
chaotic as SALI w 10~ 16 for t w 800. If Eq. flE^ is valid, the slope of the SALI in Fig. 
EJb) should be given approximately by o~i, being actually —a\j In 10, because log(SALI) 
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Figure 6. (a) The evolution of the two Lyapunov exponents an, (Jit for the chaotic 
orbit with initial condition x = 0, y = 0, z = 0, p x = 0, p y = 0, p z ~ 0.123693 
of the 3D system (JjjJ. (b) The SALI of the same orbit (solid line) and the function 
ce -(^-<T 2 )t (d as hed line) for a x = 0.0107, <r 2 = 0.0005 and c = 10 5 . Note that the 
t-axis is linear. 

is a linear function of t. As we are only interested in the slope of the SALI, we plot in 
Fig. Efb) Eq. (|T5|) for an appropriate value of c (here c = 1) and o\ = 0.047 (dashed 
line) and find that the agreement of the approximate formula (jl4j) to the computed 
values of the SALI is indeed quite satisfactory. 

Chaotic orbits of 3D Hamiltonian systems generally have two positive Lyapunov 
exponents, o\ and a 2 . So, for approximating the behavior of the SALI by Eq. ([HJ), both 
o"i and er 2 are needed. We compute <7i, a 2 for a chaotic orbit of the 3D system (0) as 
the long time estimates of some appropriate quantities an, a it by applying the method 
proposed by Benettin et al. j2J. The results are presented in Fig. Efa). The computation 
is carried out until au and a 2 t stop having high fluctuations and approach some non-zero 
values (since the orbit is chaotic), which could be considered as good approximations of 
their limits a\, a 2 - Actually for t pa 10 5 we have a\ t ~ 0.0107, a 2 t ~ 0.0005. Using these 
values as good approximations of a±, a 2 we see in Fig. Efb) that the slope of the SALI 
(solid line) is well reproduced by Eq. (|T4*|) (dashed line). Note how much more quickly 
the SALI's convergence to zero shows chaotic behavior, while the two LCEs, a± t a 2t , 
take a lot longer to reach their limit values. Moreover, the results presented in Figs. El 
and |H1 give strong evidence for the validity of Eq. (pPfj) in describing the behavior of the 
SALI for chaotic motion. So we conclude that for chaotic motion the SALI is related to 
the two largest Lyapunov exponents and decreases asymptotically as SALI oc e~^ (Tx ~ a2 ' t . 
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Figure 7. The PSS x — of the 2D Henon-Heiles system ©. The axis p y — is also 
plotted. 

4. Distinguishing between regions of order and chaos 

The SALI offers indeed an easy and efficient method for distinguishing the chaotic vs. 
ordered nature of orbits in a variety of problems. In the present section we use it for 
identifying regions of phase space where large scale ordered and chaotic motion are both 
present. 

In Fig.[7|we present a detailed plot of the x = PSS of the 2D Henon-Heiles system 
(jHJ). Regions of ordered motion, around stable periodic orbits, are seen to coexist with 
chaotic regions filled by scattered points. In order to demonstrate the effectiveness of 
the SALI method, we first consider orbits whose initial conditions lie on the line p y = 0. 
In particular we take 5000 equally spaced initial conditions on this line and compute 
the value of the SALI for each one. The results are presented in Fig. |H] where we plot 
the SALI as a function of the y coordinate of the initial condition of these orbits for 
t = 1000 (panel (a)) and t = 4000 (panel (b)). In both panels the data points are line 
connected, so that the changes of the SALI values are clearly visible. Note that there are 
intervals where the SALI has large values (e.g. larger than 10 -4 ), which correspond to 
ordered motion in the island of stability crossed by the p y = line in Fig. There also 
exist regions where the SALI has very small values (e.g. smaller than 10 -12 ) denoting 
that in these regions the motion is chaotic. These intervals correspond to the regions 
of scattered points crossed by the p y = line in Fig. Although most of the initial 
conditions give large (> 10~ 4 ) or very small (^ 10~ 12 ) values for the SALI, there also 
exist initial conditions that have intermediate values of the SALI (10 -12 < SALI ^ 10~ 4 ) 
e.g. at t = 1000 in Fig. Ufa). These initial conditions correspond to sticky chaotic orbits, 
remaining for long time intervals at the borders of islands, whose chaotic nature will be 
revealed later on. By comparing Figs. Efa) and|Hfb) it becomes evident that almost all 
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Figure 8. The values of the SALI for (a) t = 1000 and (b) t = 4000 for orbits of 
the 2D system © with initial conditions on the p y = line on the PSS (Fig. 0, as a 
function of the y coordinate of the initial condition. 

points having 1(T 12 < SALI < 1(T 4 in Fig. |^a) move downwards to very small values 
of the SALI in Fig. IHfb), while the intervals that correspond to ordered motion remain 
the same. Again in Fig. IH^b) there exist few points having intermediate values of the 
SALI, which correspond to sticky orbits whose SALI will eventually become zero. We 
note that it is not easy to define a threshold value, so that the SALI being smaller than 
this value reliably signifies chaoticity. Nevertheless, numerical experiments in several 
systems show that in general a good guess for this value could be < 1CT 4 . 

In both panels of Fig. |H1 around y ~ —0.1 there exists a group of points inside a 
big chaotic region having SALI > 10~ 4 . These points correspond to orbits with initial 
conditions inside a small stability island, which is not even visible in the PSS of Fig. 
Also the point with y = —0.2088 has very high value of the SALI (> 0.1) in both panels 
of Fig. |H1 while all its neighboring points have SALI < 10~ 9 even for t = 1000. This 
point actually corresponds to an ordered orbit inside a tiny island of stability, which 
can be revealed only after a very high magnification of this region of the PSS. So, we 
see that the systematic application of the SALI method can reveal very fine details of 
the dynamics. 

By carrying out the above analysis for points not only along a line but on the whole 
plane of the PSS, and giving to each point a color according to the value of the SALI, 
we can have a clear picture of the regions where chaotic or ordered motion occurs. The 
outcome of this procedure for the 2D Henon-Heiles system (JHJ), using a dense grid of 
initial conditions on the PSS, is presented in Fig. Ofa). The values of the logarithm of 
the SALI are divided in 4 intervals. Initial conditions having different values of the SALI 
at t = 1000 are plotted by different shades of gray: black if SALI ^ 10~ 12 , deep gray 
if 10~ 12 < SALI SC 10~ 8 , gray if 10~ 8 < SALI ^ 10~ 4 and light gray if SALI > 10" 4 . 
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Figure 9. Regions of different values of the SALI (a) on the PSS x = of the 2D 
Henon-Heiles system ©, at t = 1000 and (b) on the subspace y = 0, p y = of 
the PSS z — of the 3D system (jSJl, at t = 5000. In both frames initial conditions 
are colored black if their SALI sc: 10~ 12 , deep gray if 10~ 12 < SALI < 10~ 8 , gray if 
10~ 8 < SALI < 10~ 4 and light gray if SALI > 10~ 4 . 



Thus, in Fig. a) we clearly distinguish between light gray regions, where the motion 
is ordered and black regions, where it is chaotic. At the borders between these regions 
we find deep gray and gray points, which correspond to sticky chaotic orbits. It is 
worth-mentioning that in Fig. Ufa) we can see small islands of stability inside the large 
chaotic sea, which are not visible in the PSS of Fig. [71 like the one for y rs —0.1, p y ~ 0. 
Although Fig. Efa) was computed for only t = 1000 (like Fig. IHfa)), this time was 
sufficient for the clear revelation of small ordered regions inside the chaotic sea. 

The construction of Fig. E^a) was actually speeded up by attributing the final value 
of the SALI (at t = 1000) of an orbit to all its intersection points with the PSS, and by 
stopping the evolution of the orbit if its SALI became equal to zero for t < 1000. For a 
grid of 375 x 750 equally spaced initial conditions on the p y ^ part of the figure, we 
need about 2 hours of CPU time on a Pentium 4 2GHz PC. Although it is difficult to 
estimate, we expect that it would take considerably longer to discern the same kind of 
detail, by straightforward integration of Eqs. (jlj) for a similar grid of initial conditions. 

For 3D Hamiltonians the PSS is 4-dimensional and thus, not so useful as in the 
2D case. On the other hand, the SALI can again identify successfully regions of order 
and chaos in phase space. To see this, let us start with initial conditions on a 4d grid 
of the PSS and attribute again the final value of the SALI of an orbit to all the points 
visited by the orbit. In this way, we again find regions of order and chaos, which may 
be visualized, if we restrict our study to a subspace of the whole 4d phase space. As an 
example, we plot in Fig. Efb) the subspace y = 0, p y = of the 4d PSS z = of the 
3D system (JHJ), using the same technique as in Fig. EJa). Again we can see regions of 
ordered (colored in light gray) and chaotic motion (colored in black), as well as sticky 
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chaotic orbits (colored in deep gray and gray) at the edges of these regions. Pictures 
like the ones of Fig. El apart from presenting the regions of order and chaos, could also 
be used to estimate roughly the fraction of phase space volume occupied by chaotic or 
ordered orbits and provide good initial guesses for the location of stable periodic orbits, 
in regions where the motion is ordered. 

5. Comparison with other methods 

The results presented in the previous sections show that the SALI is a simple, efficient 
and easy to compute tool for distinguishing between ordered and chaotic motion. 
Implementing the SALI is an easy computational task as we only have to follow the 
evolution of an orbit and of two deviation vectors, computing in every time step the 
minimum norm of the difference and the addition of these vectors. In the case of 
chaotic motion the SALI eventually tends exponentially to zero, reaching rapidly very 
small values or even the limit of the accuracy of the computer. On the other hand, in 
the case of ordered motion the SALI fluctuates around non-zero values. It is exactly 
this different behavior of the SALI that makes it an ideal tool of chaos detection. The 
SALI has a clear physical meaning as zero, or a very small value of the index, signifies 
the alignment of the two deviation vectors. An advantage of the method is that the 
index ranges in a defined interval (SALI G [0, \/2}) and so very small values of the SALI 
(e. g. smaller than 10 -8 ) establish the chaotic nature of an orbit beyond any doubt. 

The SALI helps us decide the chaotic nature of orbits faster and with less 
computational effort than the estimation of the maximal LCE. This happens because 
the time needed for L t to give a clear and undoubted indication of convergence to non- 
zero values is usually much greater than the time in which the SALI becomes practically 
zero, as can be seen in Figs. IH El and El 

Many other chaos indicators have been introduced in recent years, some of which are 
compared in this section with the SALI. We also study in more detail the latest method 
we very recently became aware of, the so-called 0-1 test, introduced by Gottwald and 
Melbourne [Ej. 

The efficiency of the SALI was compared in [3] with the well-known method of the 
Fast Lyapunov Indicator (FLI) [T%1 IT^] and the method of the spectral distance D of 
spectra of stretching numbers [20J. It was shown that the SALI has comparable behavior 
to the FLI both for ordered and chaotic orbits, with the SALI being able to decide the 
nature of an orbit at least as fast as the FLI. An advantage of the SALI method with 
respect to the FLI is the fact that the SALI ranges in a given interval, with very small 
values corresponding to chaotic behavior, while the values of FLI increase in time, 
both for ordered and chaotic motion, but with different rates. So, the interpretation 
of different colors in color plots produced by the SALI method, like the ones of Fig. El 
does not depend on the integration time of the orbits, in contrast to similar plots of the 
FLI, since the range of FLI values changes as time grows. As was explained in detail in 
[Sj the computation of the SALI is much easier and faster than the computation of the 
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spectral distance D, mainly because we do not have to go through the computation of 
the spectra of stretching numbers. Also the SALI can be used to distinguish between 
order and chaos in the case of 2d maps, where the spectral distance D cannot be applied. 

Sandor et al. j2Ij comparing the Relative Lyapunov Indicator (RLI) method with 
the SALI showed that both indices have similar behaviors even in cases of weakly chaotic 
motion. The RLI is practically the absolute value of the difference of L t of two initially 
nearby orbits, with very small values of the RLI denoting ordered motion, while large 
differences between the L t 's denote chaotic behavior (for more information on the RLI 
method see [21])- We note that the computation of the RLI requires the time evolution 
of two orbits and two deviation vectors (one for each orbit), while the computation of 
the SALI is faster as we compute one orbit and two deviation vectors. 

Very recently, Gottwald and Melbourne introduced a new test for distinguishing 
ordered from chaotic behavior in deterministic dynamical systems: the 0-1 test. The 
method is quite general and can be applied directly to long time series data produced 
by the evolution of a dynamical system. In that sense, the 0-1 test is more general than 
the SALI for the computation of which we need to know the equations that govern the 
evolution of the system, as well as its variational equations. The 0-1 test uses the real 
valued function: 

p{t) = [ 0(s) cos(^(s)) ds , (16) 







where <fr(s) is, in general, any observable of the underlying dynamics and 

iP(t) = kt+ [ (p(s)ds , (17) 



o 



with k being a positive constant. By defining the mean-square-displacement of p(t): 



1 ' T 



and setting 



M(t) = hm ^ / {pit + r) - p{r)f dr , (18) 

1 — »oo ± 



logi 



one takes the limit 

K = lim K{t) , (20) 

t— >oo 

and characterizes the particular orbit as ordered if K = 0, or chaotic if K = 1. The 
justification of the 0-1 test, as well as applications of the method to some dynamical 
systems can be found in [TT| . 

In order to compare the 0-1 test with the SALI method we apply it in the cases 
of the ordered and chaotic orbits of the Henon-Heiles system (jHJ) presented in Fig. ^ 
Recall that in the case of the chaotic orbit the SALI determines the true nature of the 
orbit at t « 800 when SALI w 10~ 16 , or even at t w 400 if we consider the more loose 
condition that SALI w 10 -8 guarantees chaoticity. We consider as an observable <p(t) 
the quantity: 

<f>(t)=V(t)+Py(t), (21) 
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Figure 10. Application of the 0-1 test for the orbits of Fig. ^ K(t) oscillates with an 
amplitude decreasing towards K — for the ordered orbit (solid line), while it tends 
to K = 1 for the chaotic orbit (dashed line). 

while for the constant k we adopt the value used in jTTj, i. e. k = 1.7. The application of 
the 0-1 test requires a rather long time series of the observable <fi(t) in order to reliably 
compute firstly M(t) for T — ► oo (Eq. ITHj) and secondly K as the limit of K(t) for 
t — > oo (Eq. I20J). In our computations we set T = 90000 time units and compute K(t) 
for t G (0, 10000), which means that the particular orbit is integrated up to t = 10 5 time 
units. Although the assumption T 3> t, which is formally required for the convergence 
of K(t) (see fZj), is not fulfilled, the behavior of K(t) for the ordered and chaotic orbit 
is different, allowing us to distinguish between the two cases (Fig. [TU|) . In particular, 
K(t) increases for the chaotic orbit (dashed line in Fig. ITTH) showing a tendency to reach 
K = 1, while it tends to K = for the ordered orbit (solid line in Fig. ITUj) exhibiting 
fluctuations of slowly decreasing amplitude. From these results it is obvious that the true 
nature of the orbits is determined correctly by the 0-1 test, but the computational effort 
needed in order to be able to characterize the orbits is much higher than in the case of 
the SALI. This difference is due to, firstly, the more complicated way of estimating K, in 
comparison with the computation of the SALI, as we have to compute several integrals, 
and secondly, the fact that we must compute the particular orbit for sufficiently long 
time, in order to approximate the limits of Eqs. p8jl . PU|) . Of course, as we have already 
mentioned, the 0-1 test is more general as it can also be readily applied to time series 
data, without knowing necessarily the equations of the dynamical system. 

6. Summary 

In this paper we have applied the SALI method to distinguish between order and chaos 
in 2D and 3D autonomous Hamiltonian systems, and have also analyzed the behavior 
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of the index for chaotic orbits. Our results can be summarized as follows. 

• The SALI proves to be an ideal indicator of chaoticity independent of the dimensions 
of the system. It tends to zero for chaotic orbits, while it exhibits small fluctuations 
around non-zero values for ordered ones and so it clearly distinguishes between 
these two cases. Its advantages are its simplicity, efficiency and reliability as it can 
rapidly and accurately determine the chaotic vs ordered nature of a given orbit. 
In regions of 'stickiness', of course, along the borders of ordered motion it displays 
transient oscillations. However, once the orbit enters a large chaotic domain the 
SALI converges exponentially to zero, often at shorter times than it takes the 
maximal Lyapunov exponent to converge to its limiting value. 

• We emphasize that the main advantage of the SALI in chaotic regions is that it uses 
two deviation vectors and exploits at every step, their convergence to the unstable 
manifold from all previous steps. This allows us to show that the SALI tends to 
zero for chaotic orbits at a rate which is related to the difference of the two largest 
Lyapunov characteristic exponents o"i, 02 as SALI oc e~( ai ~ a2 ' t . By comparison, the 
computation of the maximal LCE, even though it requires only one deviation vector 
and one exponent, o~i, often takes longer to converge, since it needs to average over 
many time intervals, where the calculation of this exponent is independent from 
all previous intervals. The SALI was also proved to have similar or even better 
performance than other methods of chaos detection which were briefly discussed in 
Sec. El 

• The SALI G [0, v2] and its value characterize an orbit of being chaotic or ordered. 
Exploiting this feature of the index we have plotted detailed phase space portraits 
both for 2D and 3D Hamiltonian systems, where the chaotic and ordered regions are 
clearly distinguished. We were thus able to trace in a fast and systematic way very 
small islands of ordered motion, whose detection by traditional methods would be 
very difficult and time consuming. This approach is therefore expected to provide 
useful tools for the location of stable periodic orbits, or the computation of the phase 
space volume occupied by ordered or chaotic motion in multidimensional systems, 
where the PSS is not easily visualized, and very few other similar techniques of 
practical value are available. 
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